Stationary Velocity Distributions in Traffic Flows 
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We introduce a traffic flow model that incorporates clustering and passing. We obtain analytically 
the steady state characteristics of the flow from a Boltzmann-like equation. A single dimensionless 
parameter, R — covoto with Co the concentration, no the velocity range, and t^ 1 the passing rate, 
determines the nature of the steady state. When R <C 1, uninterrupted flow with single cars occurs. 
When R^$> 1, large clusters with average mass (m) ~ R a form, and the flux is J ~ i? -7 . The initial 
distribution of slow cars governs the statistics. When Po(v) ~ t) M as v — > 0, the scaling exponents 
are 7 = + 2), a = 1/2 when /j, > 0, and a = (fj, + + 2) when < 0. 

PACS numbers: 02.50-r, 05.40. +j, 89.40+k, 05.20.Dd 
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I. INTRODUCTION 

Traffic flows are strongly interacting many-body sys- 
tems. They also present a natural testbed for theories 
and techniques developed for physical systems such as 
kinetic theory and hydrodynamics. Traffic systems has 
been receiving much attention recently , and a number 
of approaches were suggested including fluid mechanics 
[§,D, cellular automata] I , particle hopping [jl2] p"5"[ , 
and ballistic motion |16 lq| ! Traffic networks can be 
viewed as low dimensional systems. For example, ru- 
ral traffic is intrinsically one-dimensional and urban grid 
traffic is two-dimensional. Despite this important sim- 
plifying feature, most studies in this area are numeric in 
nature. 

Ballistic models are harder to simulate than cellular 
automata and particle hopping models. However, they 
are more realistic since time and space are treated as 
continuous variables. They can also prove useful for ana- 
lytical treatment. An exactly solvable clustering process 
shows that extremal properties of the velocity distribu- 
tion determine the kinetic behavior fl(f |. However, it re- 
sults in ever-growing and ever-slowing jams with a trivial 
steady state in a finite system. In this study, we inves- 
tigate more realistic situations where fast cars can pass 
slow cars. This is motivated by and should be appli- 
cable to passing zones of one lane roadways as well as 
multilane highways. Our goal is to determine analyti- 
cally statistical properties of the flow such as the flux, 
and characterize their dependence on the intrinsic veloc- 
ity distribution. 




Fig.l Space time diagram of the traffic model. Formation of 
a cluster with two fast cars is shown to the left and formation 
of a one car cluster and its breakup due to escape is shown to 
the right. 



We start by formulating the model. Consider a one- 
dimensional traffic flow with sizeless cars ("particles") 
moving with a constant velocity. We assume that cars 
have intrinsic velocities by which they would drive on an 
empty road. Initially, cars are randomly distributed in 
space and they drive with their intrinsic velocities. How- 
ever, the presence of slower cars forces some cars to drive 
behind a slower car and therefore leads to the formation 
of clusters. Simple collision and escape mechanisms are 
implemented. When a cluster overtakes a slower clus- 
ter, a larger cluster forms. It moves with the smaller of 
the two velocities. Meanwhile, all cars in a given cluster 
may escape their respective cluster and resume driving 
with their intrinsic velocity. We assume a constant es- 
cape rate £q 1 . The actual collision and escape times are 
proportional to the car size and thus set to zero (these 
time scales should become important in heavy traffic). 

A heuristic argument suggests that a single dimension- 
less parameter underlies the steady state. Consider a 
state where the car concentration is cq, and the typical 
intrinsic velocity range is vq . Let the steady state cluster 
density be c < c , which implies the typical cluster size 
(to) = Cq/c. If large clusters form, (to) 1, then the 
overall escape rate can be estimated by (m)i^ 1 . Assum- 
ing that most collisions involve fast cars and slow clusters, 
the typical collision rate is cvq. In the steady state, the 
number of cars joining and leaving clusters should bal- 
ance and thus, crj/(cio) = v oc or c = (co/woio) 1 ^ 2 - This 
heuristic argument gives the leading behavior of the av- 
erage cluster size 



R 1 ' 2 



when R ^> 1, 



(1) 



where R is the ratio of the two elementary time scales, 
the escape time t csc = i n and the collision time t co \ — 
(co^o) -1 : 



R 



c v Q t . 



col 
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We term this dimensionless quantity the "collision num- 
ber" . For large collision numbers, large clusters occur 
according to Eq. (0), while for small collision numbers 
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the effect of collisions is small (m) = 1 + aR. Analy- 
sis of the master equations detailed below confirms this 
heuristic picture under quite general conditions. 

The rest of this paper is organized as follows. In Sec. II, 
the master equations are used to derive analytical expres- 
sions for various velocity distributions in the steady state. 
The leading behavior in the limiting cases of light and 
heavy traffic are highlighted in Sec. III. Explicit expres- 
sions are written for the special cases of uniform initial 
and final velocity distributions as well as discrete distri- 
butions in Sec. IV. We close with some open problems, a 
discussion, and possible applications. 



the final distribution, the observed quantity in real traffic 
flows, gives explicitly the intrinsic distribution. We con- 
firm that in the limit R — * oo, all clusters move with the 
minimal velocity P(v) — > 5{v), while in the limit R — > 0, 
all cars move with their intrinsic velocity P(v) — > Pq(v). 

It is convenient to transform the integral equation (|J) 
into a differential one. Consider the auxiliary function 

Q(v)=R- 1 +f dv'(v-v')P(v'), (5) 
Jo 

which gives the cluster distribution by second differenti- 
ation 



II. THEORY 

In the following, it is convenient to introduce dimen- 
sionless velocity v/vq — ► v, space xcq —>■ x, and time 
CQVot — > t variables. This rescales the escape rate t^ 1 
to the inverse collision number Let P(v,t) be the 

density of clusters moving with velocity v at time t. Ini- 
tially, isolated single cars drive with their intrinsic ve- 
locities drawn from the distribution Po(v) = P(v, t = 0). 
This intrinsic velocity distribution is normalized to unity, 
J dvPo(v) = 1. The flow is invariant under a velocity 
translation, and the minimal velocity is set to zero. 

Initially, the velocities and the positions of the particles 
are uncorrelated. Escape effectively mixes the positions 
and the velocities. Assuming that no spatial correlations 
develop, a closed master equation for the velocity distri- 
bution of clusters P(v,t) can be written 



P(v) = Q"(v). 



(6) 



dP(v,t) 
dt 



= R- 1 [P (v) -P(v,t)} 

-P(v,t) [ dv'(v-v')P(v',t). 
Jo 



(3) 



The density of slowed down cars with intrinsic velocity 
v is Po(v) — P(v,t). Such cars escape their clusters with 
rate ii -1 , and thus the escape term. Collisions occur 
with rate proportional to the velocity difference as well 
as the product of the velocity distributions. The integra- 
tion limits ensure that only collisions with slower cars are 
taken into account. 

Steady state is obtained by taking the long time limit 
t — > oo or d/dt — 0. Since we are primarily inter- 
ested in the steady state, we omit the time variable 
P(v) = P(v, t = oo). Equating the right-hand side of the 
master equation to zero, a relation between the intrin- 
sic car distribution and steady state cluster distribution 
emerges 



P(v) 



1 + R dv'(v-v')P(v') 



= Po(v). 



(4) 



Given the intrinsic velocity distribution this relation gives 
the final cluster velocity distribution only implicitly. In 
contrast, the inverse problem is simpler as knowledge of 



Thence, the steady state condition (||) reduces to the sec- 
ond order nonlinear differential equation 



Q{v)Q"(v) = R-'Poiv). 



(7) 



The boundary conditions are Q(0) = R^ 1 and Q'(Q) = 0. 
The cluster concentration is found from the cluster ve- 
locity distribution using 



dvP(v), 



and the average cluster mass is simply (m) 



(8) 



Fur- 



thermore, the average cluster velocity is obtained from 

/>OC 

(v) = c- 1 / dvvP(v). (9) 



Cars may drive with a velocity smaller than their in- 
trinsic one, and it is natural to consider the joint velocity 
distribution P(v,v'), the density of cars of intrinsic ve- 
locity v driving with velocity v' . The master equation for 
the joint distribution reads 

dP( -Q^ = - Rr x P(v, v') + (v - v')P{v)P{v') 

-P(v,v') f dv"(v' - v")P(v") (10) 



+ P(v') / dv"(v" -v')P(v,v"). 



The first term accounts for loss due to escape, while the 
rest of the terms represent changes due to collisions. For 
instance, the last term describes events where a w-car 
driving with velocity v" is further slowed down after a 
collision with a w'-cluster. One can verify that the total 
number of u-cars, 



P (v) = P{v)+ / dv'P(v,v'), (11) 
Jo 

is conserved by the evolution Eqs. (||) and (|l0|). 
At the steady state, the joint distribution satisfies 
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P(v, v')Q{v') = (v - v')P(v)P{v') + Q{v, v')P{v'), (12) 

obtained using the definition of Q(v) and the joint aux- 
iliary function 



Q(v,v') = / dw(w — v')P{v,w). 



(13) 



Although the collision number R does not appear in 
Eq. (n3) explicitly, it enters through Q(v) and P(v). 

Combining ( p"2| ) with Eqs. (|l3|), (^j, and using the re- 
lationship P(v,v') — d 2 Q(v,v')/dv' yields 



d_ 

dv 1 



,2/ t\ d Q(V,V') 



dv 1 Q{v') 



= {v- v')P{v)P{v'). (14) 



Integrating twice over v' gives the joint auxiliary function 
in terms of the single variable functions 



Q(v,v') = P(v)Q(v') 



du 



dw(v — w)P(w). (15) 



The boundary conditions Q(v,v) — -§^rQ(v, v')\ ,_ = 
were used to obtain this expression. Furthermore, in- 
tegration by parts of f dw(v — w)P(w) — f dw(v — 
w)Q"(w) gives 



Q(v,v') = P(v) 



Q(v)Q(v') 



du 



-(v- v') 



(16) 



Substituting Eq. ( |l6| ) into (|l2|) and then replacing PQ 
with R _1 Pq we find a relatively simple expression for the 
joint velocity distribution 



P(v,v') 



Po(v)P (v') 



du 



Q{v>) J v , [RQ{u)f 



(17) 



Another interesting quantity is the flux or the average 
velocity given by J — J dv[vP(v) + j Q dw wP(v,w)] . 
From the definition of the joint auxiliary function, the 
second integral is identified with Q(v,0), implying 



J= I dv[vP(v) + Q(v,0)]. 
'o 



(18) 



The integrand can be considerably simplified using 
Eq. (|l|), Q(0) = R' 1 , and Eq. (g). The term vP(v) 
cancels and we find a useful expression for the flux 



J= I dvP (v) 
>o Jo 



du 



[RQ(u)f 



(19) 



One can also ask for the actual velocity distribution of 
cars defined via 



(20) 



G(v) = P(v) + / dwP(w,v). 



Substituting the joint velocity distribution allows us to 
express the car velocity distribution via single variable 
distributions 



G(v) = P(v) 



1 + R 



dw Pq (w) 



du 



[RQ(uW 



(21) 



The car velocity distribution satisfies the normalization 
conditions 1 = JdvG(v) and J = JdvvG(v). 

In summary, for arbitrary intrinsic velocity distribu- 
tions, the entire steady state problem is reduced to the 
nonlinear second order differential equation (^). Given 
Q(v), steady state characteristics such as P(v), P(v,v'), 
J, and G(v) can be calculated using the explicit formulae 
@, ©, and ©, respectively. 



III. LIMITING CASES 

Although one cannot solve Eq. ([?]) analytically in gen- 
eral, it is still possible to obtain the leading behavior in 
the limits of R — > and R — > oo. 



A. Low Collision Numbers 

To analyze the flow characteristics in the collision- 
controlled regime, R <§; 1, we use Eq. (^) to write P(v) 
as a perturbation expansion in R: 



P(v)^P (v) 



1-R I dv'{v-v')P (v') 
o 



(22) 



In this limit, the auxiliary function is roughly constant 
RQ(v) = 1, and Eq. (17) gives the joint distribution to 
first order in R 



P(v,v')^R(v-v')P (v)P (v'). 
The final density and flux are 

c = 1 — cii?, J = Jo — Ji-R, 



(23) 



(24) 



with ci = JdvP (v) Jodv'(v - v')P (v'), J = Ah, 
Ji = M2 — Ml (M n are the moments of the intrinsic ve- 
locity distribution M n = J dvv n Po(v)). The coefficient 
Ji > equals the width of the initial velocity distribu- 
tion. This gives a simple intuitive picture: the larger the 
initial velocity fluctuations, the smaller the flux. By ei- 
ther substituting the joint velocity distribution into the 
definition of G(v), or from Eq. (21), the car velocity dis- 
tribution is 



G(v)^P (v) 



1 + R 



dv' (v — v)Pq(v') 



(25) 



As the integral is over the entire velocity range, the or- 
der R correction is positive for small v and negative for 
large v. In other words G(v) > Pq{v) when v < v c . 
The crossover velocity equals the average intrinsic veloc- 
ity v c = Jq = Mi, as seen from Eq. (p5j). 

We conclude that the collision-controlled limit is 
weakly interacting, explicit expressions for the leading 
corrections of the steady state properties are possible. 
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B. Large Collision Numbers 

The analysis in the complementary escape-controlled 
regime, R ^> 1, is more subtle since the condition 
R J Q dv'(v—v')P (v') <§; 1 is satisfied only for small veloc- 
ities. No matter how large R is, sufficiently slow cars are 
not affected by collisions, and P(v) is given by Eq. j2^ ) 
when v <C v* . The threshold velocity v* = v*(R) is 

estimated from R J™ dv(v* — v)Po(v) ~ 1. 

It is useful to consider algebraic intrinsic distributions 

P (v) = (fi + l)v» M>-1, (26) 

in the velocity range [0:1] with the prefactor ensur- 
ing unit normalization. For such distributions, the 

threshold velocity decreases with growing R accord- 
i 

ing to v* ~ R f+2 . For v ^> v*, the integral 
in Eq. (m dominates over the constant factor and 
RP(v) J Q dv'(v — v')P(v') ~ v^. Anticipating an al- 
gebraic behavior for the cluster velocity distribution, 
P(v) ~ R a v when v 3> v* , gives different answers for 
positive and negative fi. The leading behavior for »»d* 
can be summarized as follows 

( i?-i/(M+2)^-i /i < 0; 

P(v) ~ i R-iv-^lv/v*)]-* m = 0; (27) 
[ R-iv's^ 1 fi > 0. 

The small and large velocity components of P(v) match 
at the threshold velocity, P(v*) ~ Pq(v*). Careful anal- 
ysis, detailed in the following section, is needed to get 
the logarithmic corrections in the borderline case fi = 0. 
Substituting the leading asymptotic behavior of Eq. j2^ ) 
into Eq. (^), the average cluster size is found 

!^(m+i)/(^+2) M < 0; 
(i?/lni?) 1/2 M = o| (28) 
R 1 / 2 fi > 0. 

Similarly, the average cluster velocity defined in Eq. (||) 
is evaluated 

( R »/(p+2) fI<0 . 
(v)~ll/]nR fi = 0; (29) 
[ const fx > 0. 

Two distinct regimes of behavior emerge. For fi > 0, 
car-cluster collisions dominate while for fi < cluster- 
cluster collisions dominate. The scaling argument given 
in the introduction assumes the former picture, and thus 
it does not hold in general. A posteriori, one can extend 
the scaling argument to the fi < regime. The argument 
becomes involved, and we do not present it here. In- 
terestingly, in the cluster-cluster dominated regime, the 
scaling behavior for the average cluster size, (to) ~ R a 
with a = (fi + l)/(/U + 2), is identical to the kinetic scal- 
ing, (to) ~ (covot) a with the same a, found in the no 
passing limit . This suggests an analogy between the 



dimensionless collision number R = CoVoto and the di- 
mensionless time CoVot. The flux can be evaluated in a 
similar fashion using Eq. (|l9|), 

J~ v * ~ R-HT2. (30) 

Interestingly, the flux is proportional to the threshold ve- 
locity v* . As a result, the flux exponent 7 = l/(fi+ 2) 
is a regular function of fi unlike the cluster size exponent 
a. Eq. ( |30| ) is also consistent with identification of the 
crossover velocity v c with the marginal velocity v* . No 
flux reduction occurs when the intrinsic distribution is 
dominated by fast cars, i.e., in the limit fi — > 00. In the 
other extreme, the maximal flux reduction J ~ R^ 1 is 
realized when fx — ► — 1. 

The car velocity distribution is strongly enhanced in 
the low velocity limit, as seen by evaluating Eq. (^l| ) 

G(v) ~ R^v^(l -const, x v^ 1 ), v < v* . (31) 

As a check of self-consistency, one can verify that 
1 ~ f Q dvG(v), and J ~ v* ~ J" dvvG(v). Near the 
maximal velocity, the car velocity distribution ap- 
proaches the cluster distribution G(u) = -P(u)- 

In summary, as R — ► 00 the solution to the differ- 
ential equation (Q) exhibits a boundary layer structure. 
Inside the boundary layer, v < v* , the cluster velocity 
distribution is only slightly affected by collisions, while 
in the outer region v > v*, the cluster velocity distri- 
bution is much smaller than the intrinsic velocity distri- 
bution. The threshold velocity v* is determined by the 
small velocity behavior of the intrinsic velocity distri- 
bution, and for the algebraic distributions (^6|) we have 
found v* r~j R~~v+z — > 0. The behavior detailed above 
in the escape controlled limit is not restricted to purely 
algebraic distributions but is quite general. We conclude 
that a single parameter 

d 

/i = limw— lnPo(w) (32) 

v— >0 OV 

determines the behavior as R — > 00. In short, extreme 
statistics underly the escape-limited flow properties. Ad- 
ditionally, an interesting transition between a slow and a 
fast velocity dominated flow occurs at fi = 0. 



IV. EXAMPLES 

Although the above analysis is quite general, it applies 
only to the limiting values of R. To examine intermedi- 
ate behavior, it is also useful to obtain explicit solutions 
for some special cases. Below, we consider two relevant 
cases: uniform Pq(v) and P(v). We also obtain explicit 
expressions in the case of discrete velocity distributions. 
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A. Uniform Intrinsic Distribution 

We now consider the case of a uniform intrinsic distri- 
bution, Po(v) = 1 for < v < 1. This case appears to 
be the most relevant to real traffic flows since the intrin- 
sic velocity distribution should be regular near the the 
minimal velocity. Integrating QQ" = R^ 1 subject to the 
boundary conditions Q(0) = R^ 1 and Q'(0) — gives 

Q' 



\J2R 1 \n(RQ). Second integration gives 
dq 



RQ 



V2hTg 



and thus implicitly determines Q(v). 
leading behavior when R 3> 1, we find 



< TO > -\hr-5 



R_ 
ln~fi" 



1 



lniT 



(33) 

Evaluating the 
%■ (34) 



Fig. 2 shows the velocity distribution obtained numeri- 
cally using Eqs. (|f) and (|l|) for R = 10. For v < v*, 
G[v) > P (v), and for v > v* , G(v) Si P(v) < P {v). 
The calculated distributions are consistent with the pre- 
dictions, G(0) - R 1 ' 2 and v* ~ R- 1 / 2 . The car velocity 
distribution is linear near the origin in agreement with 
Eq. (§). 
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Fig. 2 Velocity distributions in the case of a uniform initial 
distribution Po(v) — 1, for R = 10. 



B. Uniform Cluster Distribution 

Consider the uniform final cluster distribution P(v) = 
c. This inverse problem is simple as all quantities can 
be obtained explicitly. From Eq. (|^), the auxiliary func- 
tion is Q{v) = R^ 1 + \cv 2 and from Eq. (Q) the initial 
distribution reads 



Po(v)=c 



1 + —Rev* 



(35) 



The overall initial concentration is unity, thereby relating 
R and c ^ 
Eq. ©, 



R and c via 1 = c + ^Rc 2 . The flux is calculated from 



■j (3 + A)VAtan- 1 v /A + A-ln(l + A) 
J = — , (36) 



3i? 



with A 



\Rc 



(3/2) 



These ex- 



y/1 + 2R/3 - 1 

plicit solutions agree with our low and high R predic- 
tions. For instance, when > 1 we find (m) ~ R 1 ^ 2 
and J ~ i?" 1 / 4 . If we look at the initial distribution, 
Pq(v) = (6/i?) 1 / 2 + 3t> 2 , then the constant part is negli- 
gible and the distribution corresponds to the [i = 2 case 
of the power-law distribution (|2^). For this case the size 
exponent is a = 1/2 and the flux exponent is 7 = 1/4, 
see ( p8| ) and (|3C]), in agreement with our findings. 

Substituting Pq(v) and Q(v) in Eq. (|l7| ) and perform- 
ing the integration gives the joint distribution 



p{vy)^2R-^ 2{v ~ v ' )[l ~ Xvvr) 



(37) 



1 + \v 12 

+2R- 1 \V 2 {1 + Xv 2 ) [tan-^A 17 ^) - tan^A 1 / 2 !/) 



A direct integration of the joint distribution confirms the 
conservation law (|Tl"|), thus providing a useful check of 
self-consistency. The joint velocity distribution is linear 
in the velocity difference for small v and v' . This is remi- 
niscent of the small collision number behavior of Eq. j2^). 
As the velocity difference increases, significant curvature 
develops (see Fig. 3). 




Fig. 3 The joint velocity distribution for the uniform final 
distribution case with R = 10. 
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C. Discrete Velocity Distribution 

The results formulated for continuous distributions can 
be used to study the special case of discrete velocity dis- 
tribution as well. Here wc quote the results in terms of 
the original (non-dimensionless) quantities. Consider the 
intrinsic velocity distribution 



Po(v) 



' CiS(v 



(38) 



with vi < V2 < ■ ■ ■ < v n . We denote by pi the dis- 
crete counterpart of the cluster velocity distribution, e.g., 
P( v ) = Y^i=\Pi^{ v ~~ v i)- The steady state condition of 
Eq. (|§) reads 



3=1 



(39) 



Substituting the intrinsic velocity distribution and solv- 
ing iteratively, we get 



Pi 

P2 
P3 = 



Cl 



C2 



1 + ci(v 2 - Vl)t 
C3 



(40) 



)'n 



etc. Rather than a solution to a differential equation, 
the steady state solution is in the form of an explicit 
continued fraction. This expression involves the initial 
distribution and the velocity differences, and can be use- 
ful to analyze data in a histogram form. In a similar 
way, explicit expressions can be obtained for the rest of 
the steady state properties. 



V. DISCUSSION 

An important property, the cluster size distribution is 
absent from our treatment so far |nj. Naturally, the 
size and the velocity of a cluster are strongly correlated 
and one must consider P m (v), the distribution of clusters 
of size to and velocity v. The joint cluster size-velocity 
distribution obeys the master equation 



dP m (v) 

at 



= R-^mPrn+^v) - (to - l)P m (v)} 
+ 5 m< i[Po(v) - P(v)} - F(v)P m (v) 

/DO m 

.7 = 1 



(41) 



which applies for all to > 1. Terms proportional to 
account for escape, while the rest represent colli- 
sions. The factor F(v) = J °° dv'\v — v'\P(v') measures 



the overall collision rate experienced by a u-cluster, and 
is reminiscent of kinetic theory. Summing Eqs. (|4l]), one 
recovers the rate equation (||) for P(v) — J2 m Pm{v). On 
the other hand, integration over the entire velocity range 
does not reduce Eqs. ([ll]) to a closed system of rate equa- 
tions for the cluster size distribution P m = J dvP m (v). 
Therefore, the entire joint distribution is needed to de- 
termine P m . Additionally, we note that in Eqs. (^|) and 
dlO|), the integration limits include only slower velocities, 
a feature that considerably simplifies the analysis. This 
property is lost for Eqs. (|4l|), thereby putting analytical 
solution out of reach. 

Nevertheless, a leading order analysis is still possible 
for low collision numbers. In the limit R <C 1, we find 



1-R / dv'\v -v'\P (v') 







P 2 (v)^RP {v) dv'(v' -v)P (v'), (42) 

J V 

P m {v) Si R^P^v). 

Heuristically, clusters with to cars are created by to — 1 
collisions and a factor R is generated in each collision. 
Although the functions P m (v) are quite complicated, the 
overall prefactor suggests an exponential cluster 

size distribution in the dilute limit. 

In the special case of a bimodal velocity distribution, 
a solution is possible. The structure of clusters here is 
simple: A cluster of size to consists of a leading slow car 
and to — 1 fast cars behind it. The rate equation ( |4l| ) 
simplifies considerably, and a Poisson size distribution is 
found P m oc e~ * f m ~ 1 / '(to — 1)!. The collision rate / is 
equal to the product of the escape time, the velocity dif- 
ference, and the fast car concentration. This steady state 
distribution satisfies a detailed balance condition as the 
escape rate and the collision rate are equal microscopi- 
cally, (to — l)P m = fPm-i- Thus, an equilibrium steady 
state is reached. However, in general, a nonequilibrium 
steady state is approached with the collision rate and the 
escape rate balancing only macroscopically. This is seen 
by noting that the cluster size may increase by an arbi- 
trary number due to collisions, but can decrease only by 
one due to escape. 

Further investigation of the collision term in the rate 
equation will be useful as well. In the no escape case 
R^ 1 = 0, the exact Boltzmann equation 



dP{v,t) 

at 



-P(v,t) I dv'(v-v')P (v') (43) 
'o 



is different from our master equation as Pq(v') replaces 
P(v', t) in the integrand [^6|. This seemingly small differ- 
ence is important as it shows that the system remembers 
the initial state. We argue that escape, no matter how 
small, induces mixing and acts to erase this memory, and 
therefore Eq. (g). It still remain, however, to establish 
quantitatively how appropriate is this mean field assump- 
tion. 



G 



The model and the results presented above can be gen- 
eralized to study other traffic situations. First, a mul- 
tilane flow can be treated as a system of coupled one 
lane flows. Escape naturally couples neighboring lanes. 
Second, a natural generalization is to heterogeneous sit- 
uations where passing is allowed only in a fraction r of 
the road. We expect that for regular distribution of these 
passing segments the problem should reduce to the homo- 
geneous case with a renormalized collision number R/r. 
The most challenging question appears to be the role 
played by the escape mechanism. We consideredthe case 
where all cars are equally likely to escape. This assump- 
tion simplified the master equation considerably as the 
escape term is linear in P(v), and thus, is exact. The 
complementary case where only the first car in the clus- 
ter can escape is interesting as well. For low collision 
numbers, large clusters are unlikely, and the behavior 
is independent of the escape mechanism. However, for 
high collision numbers the escape mechanism becomes 
weaker and larger clusters form. Indeed, a scaling ar- 
gument along the lines of Eq. ([!]) gives (m) ~ R in the 
car-cluster dominated regime. 

In conclusion, despite the simplifying assumptions 
made, the suggested model results in quite realistic be- 
havior. The overall picture is both familiar and intuitive: 
due to the presence of slower cars, clusters form and the 
overall flux is reduced. For heavy traffic, the characteris- 
tics of the flow are solely determined by the distribution 
of slow cars. A single dimensionless parameter, the col- 
lision number R ultimately determines the nature of the 
steady state. The stationary distributions obtained ana- 
lytically provide a simple practical recipe for calculating 
the flow properties for arbitrary intrinsic distributions. It 
will be interesting to analyze observed traffic data using 
these theoretical tools. 
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